Relaxation in open one-dimensional systems 
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A new master equation to mimic the dynamics of a collection of interacting random walkers in 
an open system is proposed and solved numerically. In this model, the random walkers interact 
through excluded volume interaction (single-file system); and the total number of walkers in the 
lattice can fluctuate because of exchange with a bath. In addition, the movement of the random 
walkers is biased by an external perturbation. Two models for the latter are considered: (1) an 
inverse potential (V cx where r is the distance between the center of the perturbation and the 
random walker and (2) an inverse of sixth power potential (V oc The calculated density of the 
walkers and the total energy show interesting dynamics. When the size of the system is comparable 
to the range of the perturbing field, the energy relaxation is found to be highly non-exponential. 



In this range, the system can show stretched exponential (e 



and even logarithmic time 



dependence of energy relaxation over a limited range of time. Introduction of density exchange in 
the lattice markedly weakens this non-exponentiality of the relaxation function, irrespective of the 
nature of perturbation. 



I. INTRODUCTION 

Relaxation dynamics of interacting particles in a one- 
dimensional system is often difficult to understand be- 
cause the traditional coarse-grained descriptions (such 
as hydrodynamics or time dependent mean-field type ap- 
proximations) fail in this case. This is because of the ex- 
istence of long-range correlations mediated through the 
excluded volume interaction. In such cases, random walk 
models have often proved to be successful in describ- 
ing the non-exponential relaxation commonly observed 
in one-dimensional systems. 

Recently, several experimental studies have reported 
energy relaxation in one-dimensional or quasi one- 
dimensional systems such as DNA|l|, |3(- These experi- 
ments employed the time dependent fluorescence Stokes 
shift (TDFSS) technique to gather a quantitative mea- 
sure of the time scale involved. In one of these exper- 
iments Bruns et al.\j^ have calculated structural relax- 
ation of DNA oligonucleotides. They found that the red 
shift of the fluorescent spectrum follows a logarithmic 
time dependence. The experiment described above is an 
example where dimensionality plays an important role in 
the energy and density relaxation. 

In this work, we propose a random walk model for 
the carriers in a one-dimensional channel to mimic the 
relaxation of energy and density found in experiments 
described above. There were several theoretical stud- 
ies devoted to random walk model in one-dimensional 
systems 0,0, 0,0,0,11, 13- Many recent studies based on 
random walk model in onc-dimcnsional systems have ex- 
plored transport and other dynamical properties such as, 
DC conductivity, frequency dependent conductivity, ef- 
fect of bias on diffusion, space and time dependent prob- 



ability distribution of random walkers, mean residence 
time and mean fi rst p assage time of random walkers etc. 

The work presented here is based on a master equa- 
tion for random walk in a single-file system (a one- 
dimensional system that does not allow particles to pass 
through each other) that allows exchange of the particles 
with the bath (see figure I for a schematic illustration). 
Here we employed two model potentials with different 
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FIG. 1: The perturbation potential and the various transi- 
tions are shown schematically. The attractive part of the 
Lennard- Jones potential is shown by the line (i) and Coulomb 
potential is given by the line (ii). 

characteristics to study the effect of perturbation on the 
density and energy relaxation processes. The potentials 
used here are (I) inverse of distance potential (example 
is a Coulomb's field generated by a trapped charge in 
the lattice) and (2) inverse of sixth power of distance 
potential (this is a short range interaction compared to 
Coulomb interaction and is taken as attractive part of 
the Lennard- Jones interaction). Energy and density re- 
laxation functions in a one-dimensional channel without 
number fluctuation is calculated and then compared with 
that of a channel, where particle number fluctuates due 
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to exchange with the bath. It is observed that density 
fluctuation in the channel makes energy relaxation expo- 
nential. 

The rest of this paper is organized as follows. Section 
2 gives the description of the master equation used in the 
simulation. The details of the Monte Carlo simulations 
are given in the section 3. Results of the simulations 
were analyzed in the section 4. Section 5 presents the 
conclusions from the Monte Carlo simulations. 



II. THE GENERALIZED MASTER EQUATION 

Let the total number of lattice sites in a system con- 
sisting of a linear lattice and a bath be Mr — Ml + Mb > 
where Ml is the number of lattice sites of the lattice 
and Mb is the number of accessible bath sites. The time 
dependent probability Pf(t) of finding a particle in the 
i th site of the lattice at time t is given by the following 
master equation 



The resultant master equation is simpler and is given by 



ApL(f\ J " L ' 



Ml i 



-Pt{t) £ (t) + £ Wj 4t)pf(t), (i) 

3=1 3=1 

(the prime on the summation signifies that the term 
i = m is to be omitted from the sum) where w m ^(t) gives 
the transition probability of the particle from site m to i 
per unit time. The last two terms of the master equation 
introduce density fluctuation in the lattice by allowing 
exchange of particles with the bath. In the master equa- 
tion, the index j sums over the bath sites and P B j is 
the probability of finding a particle in the j th bath site. 
This master equation can be simplified by assuming that 
lattice sites in the bath are numerous ( Mb > > Ml ) and 
the number of walkers in the bath is much larger than 
that in the lattice. These assumptions allow us to per- 
form an averaging over the bath sites. Summation in the 
last two terms of the equation (2) can be omitted by sub- 
stituting the properties of the bath with that of an ideal 
bath defined below. The sum of the transition probabil- 
ity Wi.j(t) to the bath sites can be replaced by the term 
= w out (t)), since Wij(t) depends 
only on the temperature, size and number of particles 
in the system (discussed in detail later in this section). 
This ideal bath has infinite capacity to absorb particles. 
Similarly, when a particle is absorbed by the lattice this 
bath site acts as a supplier of infinite number of parti- 
cles. The transition probabilities from all lattice sites to 
the bath are equal and assumed to be independent of the 
instantaneous state of the bath. Therefore, Wjj can be 



replaced by w in and Y^=i Pf can be replaced by Pi 



d_pm 

dt 



= W m At)P^)-^n{t)Pt(t) 
m 

-WoutPi{t) + W in {t)P bat h- 



(2) 



If only nearest neighbor exchanges are allowed then the 
master equation can be written as 



dpm 

dt 



{t)Pt+i (*) + Wi-l,i {t)Pt-i (*) 

-Wi ti+ i(t)Pf(t) - W i ,i-i(t)Pf , (jt) 

-W out {t)Pt{t)+w m {t)P bath . (3) 



At equilibrium this system represents a grand canonical 
ensemble. Hence for insertion and removal of the parti- 
cles in the randomly selected sites we use grand-canonical 
Monte Carlo method, where probability (P u ,l,t) of sys- 
tem having number of particles N is|18j 



L 



N 



A N Nl 



-P[E-nN] 



(4) 



where j3 is the inverse of Boltzmann constant times ab- 
solute temperature {j^f), A (= ^/ 2 7rmfc B f ) ^ s ^ e ther- 
mal de-Broglie wavelength, L is the length of the linear 
lattice, fi is chemical potential and E is the potential en- 
ergy of the linear lattice. Particles of this system obey 
Boltzmann distribution at equilibrium, hence the tran- 
sition probability u^i+i (t) can be calculated from the 
total energy cost for the move. The hopping probability 
is calculated for a particle from one site to another using 
Metropolis scheme 01 • Hopping probability of a random 
walker from one site to the neighboring site of the lattice 
is given by 



u i,i+ 



i = mffl[l, e- 0AE \. 



(5) 



The transition probabilities for absorption and desorp- 
tion from the lattice can be obtained [18| so as to satisfy 
the detailed balance at equilibrium as 



Wi n — min[l, 



Q 



(iV + l) 



-P AE] 



N 



Wout = min[l, —e 13 ], 



(6) 
(7) 



bath ■ 



where Q is given by £ e^. A non-trivial aspect of this 
master equation is that the transition probability Wi n (t) 
varies with time in response to the number fluctuations. 
This explicit time dependence of Wi n (t) pose formidable 
difficulty in obtaining an explicit analytic solution of the 
problem. In this system of interacting particles, the tran- 
sition probability to neighboring sites depends on the in- 
stantaneous probability of that site being occupied. In a 
one-dimensional channel with hard rod interactions be- 
tween the carriers, mobility of each carrier is restricted to 
a portion of the lattice. The particles move in the linear 
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lattice under the influence of a chosen potential, which is 
at a fixed position of this lattice. Hence the energy cost 
for hopping in the lattice is 



AE 



Ei (transition to vacant site) 
oo (transition to occupied site) 



where E± = Ki(-^ — ^— -) for the potential 1 and E\ = 
K 2 {-^ — a ,. 1 i ( i ) for the potential 2. For absorption of 
particles to the lattice at ith site 



AE 



E 2 (creation of vacant site) 
oo (creation of occupied site) 



where E 2 — K\±- for potential 1 and E 2 = K 2 -^ for 
potential 2. For the desorption of particles from the ith 
site 



AE = E 2 



(8) 



Q in this system is constant and equal to the initial num- 
ber of the particles in the lattice. This assumption spec- 
ifies the value of the chemical potential of this linear lat- 
tice. 



constant K 2 is kept high to extend the range of the po- 
tential, value of K 2 used is -10 k B T. 

The non-equilibrium Monte Carlo simulation starts 
from a randomly chosen initial configuration. Then a per- 
turbing potential is introduced in the lattice at at time 
t = 0. Subsequent relaxation of energy is monitored. 
The simulation is repeated with different initial configu- 
rations and the results are averaged. In the Monte Carlo 
simulation movements are chosen randomly such that no 
two events can occur at the same time. In this simula- 
tion one Monte Carlo step is equivalent to one unit of 
time. Total potential energy of the lattice at any instant 
of time when perturbed by Coulomb potential is given as 



S(0 = ^M-Pi(t), 



(9) 



where K\ is a constant. Similarly, for the simulations that 
use attractive Lennard-Jones potential as the perturba- 
tion, the instantaneous potential energy of the lattice is 
given by 



(10) 



III. THE DETAILS OF MONTE CARLO 
SIMULATION 

The Monte Carlo simulations are carried out on a lin- 
ear lattice, with 50% (on the average) of the lattice sites 
are occupied by the random walkers. Closed boundary 
condition is used to explore the size effect on the relax- 
ation. In this one-dimensional lattice, the adjacent sites 
are equally spaced and all the sites are identical in energy 
in the absence of the perturbation. Here inhomogencity 
in the lattice is generated by the perturbation of potential 
introduced in the lattice. In addition, there is a dynamic 
disorder in the lattice that originates from the instanta- 
neous rearrangement of the interacting random walkers. 
The initial configuration of the random walkers in an un- 
perturbed lattice is chosen from a random distribution, 
such that no two particles occupy the same site. The in- 
verse potential arise from Coulomb interaction, hence the 
magnitude of the biasing potential is calculated from the 
interaction between the charge of the carriers (qi) and 
the charge at the center of bias (q 2 ) in a medium of di- 
electric constant e. The constant K\ for potential energy 
can be calculated as K\ = 3iM_ Here both the perturb- 
ing charge and the charge of the carrier have magnitude 
of one electron and their signs are opposite. Assuming 
high screening effect, the value of e is taken as equal 
to that of water, that is 80. The distance between the 
adjacent sites in the linear lattice is 4A(approximately 
the vertical distance between two DNA base pairs in the 
double helix). All simulations were conducted near room 
temperature (300K). The value of K x used in the simu- 
lation is -l.7ksT. Compared to Coulomb's interaction, 
the Lennard-Jones interaction is short ranged, hence the 



where k 2 is a constant. 

For comparison of relaxation function for different sys- 
tem sizes, we calculate dimensionless energy relaxation 
function S(t) which can be defined as 



S(t) = 



E(t)-E(oo) 
E(0) - E(oo) : 



(11) 



where E(t) is the instantaneous energy at time t and 
E(oo) is the average energy of the lattice, in equilib- 
rium with perturbation. Density relaxation is measured 
in terms of a dimensionless quantity 



C(t) = 



X(t)-X(oo) 
X(0)-X(oo)' 



(12) 



where X(t) is defined as X(t) = J2i XiPi{t). 

In the equilibrium simulation the system is allowed to 
equilibrate with perturbation for 10,000 steps to get the 
initial configuration. The fluctuations in total energy of 
the system can be defined as F(t) — < E(0)E(t) >. 
Here a dimensionless energy fluctuation relaxation func- 
tion can be defined as 



S(t) 



F(t) - F(oo) 
F(0) -F(oo)' 



(13) 



IV. RESULTS AND DISCUSSIONS 

A. Relaxation under Coulomb potential 

Figure 2 shows the energy relaxation function for dif- 
ferent system sizes obtained from the non-equilibrium 
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FIG. 2: Energy relaxation function S(t) is plotted for four dif- 
ferent system sizes under Coulomb potential: Ml =50 (dashed 
line), 7Vl=100 (dotted line), 7Vl=150 (dash-dot line) and 
A/i=200 (continuous line)in a closed system. The inset shows 
the fit of the energy relaxation function for A/l=200 using 
function fi(t) (dotted line) and f2(t) (dash-dot line). 
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FIG. 4: Energy relaxation function is plotted for different 
system lengths (A/l=50 (continuous line), JVz, = 100 (dash-dot 
line), A/"l=150 (dotted line) and A/l=200 (dashed line)) under 
the perturbation of Coulomb potential in an open system. 
The log of S(t) versus t plot in the inset shows straight lines 
due to the exponential relaxation. 
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FIG. 3: The density relaxation function C(t) is plotted in a 
closed system for four system sizes under Coulomb potential: 
7Vl=50 (dashed line), JVl = W0 (dotted line), A/"l=150 (dash- 
dot line) and A/i=200 (continuous line). The inset shows 
double exponential fit (continuous line) for C(t) at A/l=200. 



simulations of a one-dimensional channel without par- 
ticle fluctuation. At short times, the potential energy 
of the lattice relaxes faster, as a response to the newly 
created center of perturbation. This results in the accu- 
mulation of random walkers near the center of the biasing 
field, which slows down the relaxation rate. At the same 
time, channels far from the center of perturbation remain 
active due to the reduction in the density of particles in 



500 1000 1500 2000 

time 

FIG. 5: Density relaxation function with density fluctuation 
C(i) is plotted here for four system sizes (A/l=50 (continu- 
ous line), 7Vl=100 (dash-dot line), A/i=150 (dotted line) and 
A/l=200 (dashed line)) under Coulomb potential. 



that region. These effects together contribute to the ob- 
served slowing down of the relaxation process after the 
fast initial decay. As the system size increases the lowest 
possible energy accessible to the system obviously be- 
come lower, hence the relaxation becomes progressively 
slow as the system size increases. However, this effect be- 
comes insignificant beyond a limiting size of the lattice 
beyond which the strength of perturbation becomes neg- 
ligible. To analyze the behavior of relaxation functions 
in the system, we have fitted the relaxation function to 
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FIG. 6: Energy fluctuation relaxation function is plotted 
for a closed system at four sizes (A/"l=50 (continuous line), 
yV L = 100(dash-dot line), A/i=150 (dotted line) and 7Vl=200 
(dashed line)) in equilibrium under the perturbation of 
Coulomb potential. The inset shows the fit of the energy 
relaxation function for A/l=200 using function fi(t) (dotted 
line) and ,f2(t) (dash-dot line). 

a sum of an exponential and a stretched exponential of 
the form 

h{t) = b 1 e- t ^ +b 2 e-^l T ^ (14) 

(with constraints 61+62 = 1 and < 61, 62 < 1). In order 
to compare with the nature of relaxation with the results 
obtained by Bruns et al. 0] , a logarithmic function of the 
form 

h{t) = 1 - ci log w (at) + c 2 e-*^ (15) 

is also used to fit the relaxation function. The inset of 
figure 2 shows the fit of the energy relaxation function 
with fi(t) and / 2 (i), when size of the system is A/"l=200. 
Fitting parameters obtained for the function fi(t) are 
h = 0.20, r x = 1.4 x 10 5 , r 2 = 3.4 x 10 5 , (3 = 0.52. The 
fitting parameters for the function f%{t) are c\ = 0.16 
a = 2.5 c 2 = 0.25 and 7 = 7.1 x 10 3 . It is evident 
from the figure that the stretched exponential function 
and logarithmic functions both give good description of 
the energy relaxation function in a one dimensional lat- 
tice without density fluctuations, in the presence of the 
Coulomb potential. The time dependence of stretched 
exponential |l7l Il9| relaxation, is given by the function 
of the form S(t) = S e-^ T "^ , where < /3 < 1. Theo- 
retical explanation of the origin of stretched exponential 
relaxation in the condensed matter have been addressed 
by many |20l l2ll 1221,1231, including Huber and coworkers 
in a series of papers |2^[25| . Their model is based on the 
following simple master equation approach 

= w ^P m {t) - w itfn Pi(t) (16) 



FIG. 7: Dimensionless energy relaxation function S(t) of a 
one- dimensional channel with constant number density at dif- 
ferent system sizes (A/"l=50 (continuous line), A/l=100 (dash- 
dot line), ./Vl = 150 (dotted line) and _/V L =200 (dashed line)) 
under the perturbation of Lennard- Jones potential is shown 
here. The inset shows the fit of the energy relaxation func- 
tion for A/l=200 using function fi(t) (dotted line) and f2(t) 
(dash-dot line). 
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FIG. 8: Density relaxation function of the lattice under 
the perturbation of Lennard- Jones potential at different sys- 
tem sizes (A/l=50 (continuous line), A/l=100 (dash-dot line), 
A/l = 150 (dotted line) and A/l=200 (dashed line)) for a closed 
system is shown in the figure. 



where w m .i gives the transition probability of the particle 
from site m to i. Stretched exponential relaxation can 
arise when there is a continuum of relaxation channels 
and the probability of any single channel being open is 
much less than unity. 

Figure 3 shows the density relaxation function (which 
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FIG. 9: Energy relaxation function (S(t)) is plotted at 
different system sizes (Al=50 (continuous line), A/l = 100 
(dash-dot line), 7Vl=150 (dotted line) and A/l=200 (dashed 
line)) with particle number fluctuations under Lennard- Jones 
potential. The log of S(t) versus t plot in the inset shows 
straight lines due to the exponential relaxation. 
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FIG. fO: Energy fluctuation relaxation function is plotted 
for system in equilibrium under the perturbation of Lennard- 
Jones potential for different system sizes (A/l=50 (continu- 
ous line), A/*l = 100 (dash-dot line), A/l=150 (dotted line) and 
A/l=200 (dashed line)). The inset shows the fit of the energy 
relaxation function for A/l=200 using function fi(t) (dotted 
line) and f2(t) (dash-dot line) 



is a measure of particle aggregation) for different sizes 
of a system without number density fluctuations. At 
A/l=200, a double exponential 



h(t) = cr e-'^ 1 + c 2 e 



-t/72 



(17) 



(with constraints c\ + c 2 = 1 and < ci, c 2 < 1) (shown 
in the inset of figure 3) fit of the density relaxation func- 



tion reveals the presence of two distinct time scales in 
the density relaxation. The fitting parameters obtained 
for the double exponential { fzit)) fit of density relaxation 
function arc c x = 0.90, 71 = 1.4 x 10 6 , 72 = 8.3 x 10 4 . 

We now turn to the case when the number of walkers 
in the lattice can fluctuate due to the exchange with the 
bath. Figure 4 gives the energy relaxation function plot- 
ted for different system sizes for this case. Note that, in 
this case, the energy relaxation is faster than the non- 
fluctuating case. The exchange of particles with bath is 
equivalent to opening up of many wider channels of re- 
laxation that dominate the relaxation process. The log 
S(t) versus t plot in the inset of figure 4 show straight 
lines, which is an evidence of exponential relaxation. In 
this case walkers can bypass the obstacle caused by hard 
rod interaction in the path by exchange of particles with 
the bath. Figure 5 shows the density relaxation function 
of an open system for different sizes. As system size in- 
creases, the effect of Coulomb field on the distribution 
of carriers decreases. Hence the carrier density oscillates 
around a mean value due to the density fluctuations and 
the effect of perturbation in the density relaxation re- 
mains feeble and short lived. 

The energy fluctuation relaxation function for differ- 
ent sizes of a system with constant number density and 
that is in equilibrium with perturbation is shown in the 
figure 6. This relaxation function shows very slow non- 
exponential decay. Inset of figure 6 shows the energy 
relaxation function obtained is well fitted by f\(t) and 
f 2 (t). The fitting parameters obtained for the /i(i) are 
hi = 0.2, n = 3.7 x 10 3 , r 2 = 3.6 x 10 4 , /3 = 0.34 and 
the fitting parameters obtained for f 2 (t) are ci = 0.21 
a = 0.06 c 2 = 0.13 and 7 = 6.9 x 10 3 . Note that the en- 
ergy relaxation function in non-equilibrium and energy 
fluctuation relaxation function in equilibrium show loga- 
rithmic time dependence. 



B. Relaxation under Lennard- Jones potential 

The pronounced non-exponential decay found in non- 
equilibrium energy relaxation function for Coulomb po- 
tential is also found in the case of attractive Lennard- 
Jones potential. Figure 7 shows the energy relaxation 
function at different sizes for a system without num- 
ber fluctuation perturbed by the attractive Lennard- 
Jones potential. The fit of the energy relaxation func- 
tion with fi(t) and f 2 (t) are shown in the inset of the 
figure 7. The fitting parameters obtained for fi(t) arc 
61 ~ 0, 6 2 ~ 1.0, t s = 1.1 x 10 3 , (3 = 0.62 and f 2 (t) are 
ci = 0.14 a = 5.6 x 10 1 c 2 = 0.58 and 7 = 1.3 x 10 3 . It is 
evident from the figure that, in the case of short-ranged 
interaction, the appropriate function which can fit re- 
laxation function is stretched exponential. Here due to 
the large strength of the potential near the perturbation 
center, the initial relaxation is driven and faster. Corre- 
sponding density relaxation function for different sizes of 
this system is shown in the figure 8. 
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The energy relaxation function (S(t) versus t) of a lin- 
ear lattice with particle number fluctuation is plotted for 
different system sizes in the figure 9. The inset of the 
figure 9 shows log of S(t) versus t plot of the energy 
relaxation function, which show straight lines that is a 
signature of exponential relaxation. It is clear from the 
figures 4 and 9 that the number fluctuation in the system 
makes the energy relaxation nearly exponential irrespec- 
tive of the nature and range of the perturbing potential. 

The equilibrium energy fluctuation relaxation function 
of a system under the Lennard-Jones perturbation, for 
different system sizes, is plotted in the figure 10. The 
energy fluctuation in this system is mostly from the re- 
gion where potential is weak. The fit of the energy 
fluctuation relaxation with functions fi(t) and f2(t) are 
shown in the inset of the figure 10. The fitting pa- 
rameters for the function f\(t) are b\ = 0.42, t\ = 
8.8 x 10 3 , t 2 = 1.6 x 10 5 , (3 = 0.73 and for f 2 (t) are 
Cl = 0.14 a = 5.6 x 10 1 c 2 = 0.58 and 7 = 1.3 x 10 3 . 

Finally, note that the stretched exponential fit of relax- 
ation function is more appropriate in this case due to the 
short range of the potential. We have found no evidence 
of logarithmic time dependence for the energy relaxation 
in this case. 



V. CONCLUSION 

Let us first summarize the main results of this work. 
We have demonstrated that the cooperative dynamics 
of random walkers in a simple one-dimensional channel 
can give rise to highly non-exponential relaxation, when 
number density fluctuations are not allowed. In these 
simulations, two perturbing potentials (the Coulomb and 
the Lennard-Jones) have been used to study the effects 
of perturbing potential on the relaxation process. The 
energy relaxation under these potentials can be approxi- 
mately described by a stretched exponential (in general) 
in a closed system. The variation in the time scale of 



relaxation under these two well-known potentials can be 
understood in terms of the difference in the range of these 
potentials. The Coulomb potential being long ranged (in 
comparison with the Lennard-Jones interaction), shows 
much stronger non-exponcntiality in the energy relax- 
ation. Under the Coulomb potential, the exponent (3 of 
non-equilibrium energy relaxation function is 0.52 while 
under Lennard-Jones potential it is 0.73. The simula- 
tions seem to agree with the results of Bruns et al. Q in 
showing a logarithmic time dependence of the energy re- 
laxation function under the Coulomb potential. However, 
under the short-ranged Lennard Jones potential, the en- 
ergy relaxation function does not show logarithmic time 
dependence. 

In the smaller sized systems the relaxation is found 
to be faster and as the system size increases, relaxation 
slows down, as expected. Also as expected, the density 
fluctuations in this one-dimensional channel make the re- 
laxation function faster and exponential. When number 
fluctuation is allowed, random walkers overcome the re- 
sistance of the hard rod interactions, by moving in and 
out of the linear lattice, such that the random walkers 
experience no major hindrance to their flow. In this case 
the system behaves like a system of weakly interacting 
particles, without the need for strong cooperativity for 
relaxation. It is worth noting that interactions found in 
nature are numerous, but the models of asymptotic time 
dependence of relaxation commonly found in nature are 
limited in number. However, the parameters of the relax- 
ation function depend on the form of the biasing potential 
and the nature of interaction between the carriers. 
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